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Abstract 

We perform Monte-Carlo simulations of the three-dimensional Ising model at the 
critical temperature and zero magnetic held. We simulate the system in a ball with 
free boundary conditions on the two dimensional spherical boundary. Our results for 
one and two point functions in this geometry are consistent with the predictions from 
the conjectured conformal symmetry of the critical Ising model. 
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1 Introduction 

It is a long standing conjecture that continuous phase transitions are described by conformal 
invariant field theories [T]. Under reasonable assumptions this conjecture has been proven 
in two dimensions PE] and, recently, in four dimensions PE]. A general proof in three 
dimensions has not yet been found. However, assuming the validity of this conjecture it has 
been possible to formulate conformal bootstrap equations and find approximate solutions 
that predict the Ising critical exponents with high accuracy PEI El El no] . This success 
can also be viewed as strong evidence for conformal invariance of the 3D Ising model at the 
critical temperature. 

Conformal invariance of the critical 3D Ising model can also be tested directly with lattice 
Monte-Carlo simulations. In particular, in this work, we test the predictions of conformal 
symmetry for the critical Ising model in a ball with free boundary conditions on the two 
dimensional spherical boundary. In two dimensions, the analogous geometry (disk) was an¬ 
alyzed in mi. There have been other Monte-Carlo studies of conformal invariance in the 
3D Ising model. Using the standard cubic lattice hamiltonian, ra showed that some two 
point functions in the presence of a line defect have the functional form predicted by con¬ 
formal invariance. In p3], the authors used an anisotropic hamiltonian (with a continuous 
direction) to simulate the 3D Ising model in several cylindrical geometries and measured cor¬ 
relation functions compatible with conformal invariance. See also mm for an alternative 
implementation of a 3D cylindrical geometry. 


2 Ising model and Conformal Field Theory 

The Ising hamiltonian is 

#[{»}] = - £ s ( x ) s (y) (i) 

<x,y> 
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where the local spin variables can take two values s(x) = ±1 and the sum is over nearest 
neighbours in a cubic lattice. We are interested in correlation functions of local operators at 
the critical temperature 


{Ol{xi) . . .On{x n )) = 4y^ e ^ [{s}] Oi(zi)...O n (£ n ) , (2) 

z (4 

where the partition function is 

Z = e~ pcH[{s}] . (3) 

(4 

The local operators in the Ising model can be classified by their quantum numbers with 
respect to the Z 2 spin-flip symmetry and the point group symmetry of the cubic lattice. The 
simplest local lattice operator that is invariant under the lattice symmetries that preserve 
the point x and is odd under spin flip is the local spin held s(x). In the sector of operators 
invariant under the spin-hip symmetry and the lattice symmetries the simplest local operators 
are the identity / and the energy density 

e (x) = ^ ^2 s{x)s(x + 5) (4) 
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where x + 6 runs over the 6 nearest neighbours of x. 

At the critical temperature, the Ising model has infinite correlation length and its corre¬ 
lation functions decay as power laws of the distances \xi — Xj\. The conjecture we want to 
test is that we can dehne a Conformal Field Theory (CFT) that describes the correlators for 
|Xi — Xj much greater than the lattice spacing a. In this continuum limit, the local lattice 
operators can be written in terms of the operators of the Ising CFT that have the same 
symmetry properties. In particular, the spin held can be expanded in terms of the CFT 
scalar operators a, a', .. .which are odd under spin-hip Q 

s(T) = b sa a Aa cr(x) + b s<T ’ a A °'a'(x) + ... , A^ < A CT / < ... (5) 

where Aq is the scaling dimension of the operator O and the 6’s are dimensionless constants 
that depend on the normalization convention of the CFT operators. Similarly, the local 
energy density operator can be written as 

e(x) — b e j I + b ee a Ae e(x) + b ee > a At 'e'(x) + ... , A e < A e / < ... (6) 

where I is the identity operator and e, 4,. • • are H ie lowest dimension scalar primary oper¬ 
ators in the Z 2 -even sector. The best estimates for these scaling dimensions are 


A a = 0.518151(6) , A a r > 4.5 , A e = 1.41264(6) , A e , = 3.8303(18) . (7) 


4n general this expansion also includes descendant scalar operators, which we did not write to avoid 
cluttering. In particular, the operator d 2 a is also present and has lower dimension than o' . 
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These have been determined by a variety of methods^] like direct experimental measurements, 
Monte-Carlo simulation, high-temperature expansions, e-expansion and, more recently, con¬ 
formal bootstrap techniques. 

We normalize the CFT operators imposing the following correlation functions in the 
infinite system without boundaries 

(Oi(x)0,(y)) = , (0(x))= 0, VO#/. (8) 

This can be used to fix the coefficients in (J5| and (g- The recent Monte Carlo simulations 
of |T7] found 

b sa = 0.550(4) , b eI = 0.330213(12) , b ee = 0.237(3) . (9) 


3 CFT inside a sphere 


Let us consider the three dimensional critical Ising model in a ball with free boundary 
conditions on the two dimensional spherical boundary. We would like to test if, in the 
continuum limit, this system is described by a Boundary Conformal Field Theory (BCFT). 
In order to do this, we start by discussing what are the implications of conformal symmetry 
for correlation functions in this geometry. 


Consider first the one-point function ( 0(r )) of a scalar operator placed at distance r from 
the centre of a sphere of radius R. Without conformal invariance, this one point function 
would take the general form 



( 10 ) 


for an arbitrary function /. Imposing conformal symmetry leads to 


/(*) 


do 

(1 ^z 2 ) A ’ 


( 11 ) 


for some constant ao- 


In order to see how this comes about, we define the ball 


(x 1 ) 2 + (x 2 ) 2 + (x 3 * * * * ) 2 < R 2 , 
and introduce new coordinates y 11 via [^] 

2 R 2 y M + y 2, R 


x 11 = (R — r 


R A + 2 R 2 y ■ r + r 2 y 2 


e M 3 


+ R 


( 12 ) 


(13) 


2 See [12 for a review. 

3 It might be helpful to understand this change of coordinates as a sequence of simpler steps. Start by 

doing a translation to bring the point x 11 = to the origin and then perform an inversion. More precisely, 

defining w 11 = (R 2 — r 2 ){x M — rR/(x^ — r^) 2 maps the ball x 2 < R 2 to the exterior of a sphere of radius R 

centred at u/ M = r M . Moreover, the point x M is mapped to = oo. Then, we just translate the sphere 

in the w coordinates to the origin and perform another inversion, w M = R?y^ jy 2 . 
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where is a constant vector with norm smaller than R. This coordinate transformation 
has several nice properties. The first is that it preserves the spherical boundary. In other 
words, the ball (12) corresponds exactly to the same region in the //-coordinates 

(y 1 ) 2 + ( y 2 ) 2 + (y 3 ) 2 < R 2 , ^ e M 3 . (14) 


The second nice property is that the origin in the //-coordinates is mapped to the point 
x fj. _ inside the sphere, which we are free to choose. The third nice property is that the 
original flat metric becomes 

d4/ td2 _ 2\2 

ds 2 = d/jRdXn = ——-—- dy^dya = f 'l 2 (y)dy fl dy Pj , (15) 

M (R 4 + 2 R 2 y ■ r + r 2 y 2 ) 2 y v ’ 

in the //-coordinates, i.e. it is a conformal transformation. 

Correlation functions inside a sphere with a flat metric are equal to correlation functions 
inside the same sphere but with metric ds 2 = f l?(y)dy^dy^. 


. . . (D n (x n ))dx^dx M (^l(?/l) • • • On{]Jri)')Q. 2 (y)dyi 1 dy l j, ■ (16) 

Notice that this is true in any theory because it follows just from a relabelling of points 
without changing the physical geometry. Remarkably, correlation functions of scalar primary 
operators in CFTs also satisfy 

(0\{yi )... O n (y n )}^( y ) dy ^ dyil = Ai (l/i) • • • H An (yn){Oi(yi)... O n {y n ))d y ^d y „ • (17) 


In other words, CFT correlation functions transform in a simple way under Weyl trans¬ 
formations (or local rescalings) of the metric. Equations (16) and Jl7|) together lead to [^] 


(CM*!) • • • OnM) = FT Al (//i) • • • n~^(y n )(0 1 ( 7 / 1 )... O n (y n )) , 




for any parameter in (13). In the case of the one-point function, one finds 

A 


(0(x = r))=Q- A (0)(0(y = 0)) = 


R 1 


cio 
~R^ ’ 


R 2 — r 2 

as anticipated in ([TTj). 

For a connected two-point function inside the sphere, choosing r — x 1 , we obtain 

1—A 


(18) 


(19) 


(0( Xl )0(x 2 )) c = (0(x 1 )0(x 2 )) - (0( Xl ))(0(x 2 )) = [fi(0 )n(y)]~* (O(0)O(y)) c (20) 


where 




x 


\) x 2 — {R 2 — 2 xi ■ x 2 + xl)xi 

r.2~2 


R 4 — 2R 2 xi ■ x 2 + x \ x 2 ' 

In this way we relate a generic two point function inside the sphere to a two point function 
where one of the points is at the centre of the sphere. From spherical and scaling symmetry, 
it follows that 

(O[0)O[y)) c = ,„/ (C) ,' A . C = , (22) 


C R 2 - y 2 ) 


y 


4 From now on, we drop the subscript indicating the metric when it is the standard flat Cartesian metric. 
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Therefore, we conclude that the two point function of scalar primary operators inside a 
sphere is given by 


where 


(0( Xl )0(x 2 )) c 


r 2A 

{R 2 -x 2 1 ) A (R 2 


~TjxfooiO 


(R 2 - xj)(R 2 - x 2 ) 
R 2 (x i - x 2 ) 2 


(23) 


(24) 


The (( —y oo limit corresponds to the two points approaching each other (xi —)■ x 2 ) and 
it is controlled by the same singularity as the two point function ([8]) of the infinite system. 
This gives 

foo{C) ~ C A , C^oo. (25) 

The () —^ 0 limit corresponds to one point approaching the spherical boundary. This limit is 
controlled the boundary Operator Product Expansion (OPE) [T8I [12) 20. [2T1J 


0(z,z) 


0-0 


+ 


l oo 


(2 ;) a (2c) 


A-A 


o(z) + ... 


(26) 


where we considered a flat boundary, and used coordinates z along the boundary and the 
distance to the boundary z. Normalizing the boundary operators O to have unit two point 
function, {0(z\)0{z 2 )) = \z\ — z 2 | _2A , we obtain 


foo(C) 



C “^ 0 > 


(27) 


where O is the boundary operator with lowest dimension that appears in the boundary OPE 
of O (excluding the identity). 

We will consider the Ising model with free boundary conditions which is known to be 
described by a BCFT usually called the ordinary transition [T 8 l 1122 ]. This BCFT can be 
defined by the property that it only has one relevant boundary operator a. This operator is 
Z 2 odd and has 1.27 [231 ED 23 125]- The recent conformal bootstrap study j2B] also 
obtained El 

= 1.276(2) , a 2 ~ = 0.755(13) , a e = -0.751(4) . (28) 

The lowest dimension boundary operator in the Z 2 even sector is the displacement operator 
EH- The displacement operator D has protected dimension A^ = 3 and its correlation 
functions obey Ward identities. In particular, 

J dw(D(w)0(z,z )) = 0(z,z )) • (29) 

The two point function of a bulk and a boundary operator is fixed by conformal symmetry. 
In particular, 

(D(w)e(z, z)) = a e5 C B - —3 , (30) 

[z 2 + (z — W) z \ 

5 To estimate a e we used the value of the bulk OPE coefficient \ aae reported in [S] . 
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where Cg is the normalization of the two-point function of the displacement operator 


(D(w)D(z)) = 


C 


D 


Z — W 


i3 ' 


( 31 ) 


Using (|30|) and the Ward identity (|29|), we conclude that 47 t a e gC]g = A e a e . Finally, we 

1 ° 


conclude that 


UC)'*C 5 a' e 3 C i = 7 r; 


A f a f 


C b \ 4tt 




C^o . 


(32) 


4 Results from Monte-Carlo simulation 


In order to perform a Monte-Carlo simulation of the critical Ising model, we need to know 
the critical temperature with high precision. [^] We used Wolff’s cluster algorithm [28] to 
avoid critical slowing down and used the value 


P c = 0.22165463(8) . (33) 

of the critical temperature from [29]. To check that this is a good estimate of the critical 
temperature we measured the Binder cumulant 

Ub = 1~ ~ 2 , (34) 

3(m 2 ) 2 

where m = ~^Yl x s (x) is the magnetization per spin, with N = (L/a) 3 the total number 
of spins in a system with periodic boundary conditions. In figure [lj we plot the Binder 
cumulant for several system sizes. 


We also used these simulations with periodic boundary conditions to determine the ex¬ 
pectation value ( e(x )) = b e j of the energy density operator in the infinite system at j3 = /3 C . 
In figure [2j we show the measured values of {e{x))i for several system sizes and fit them to 
the theoretical expectation from equation 


{e(x)) L = b eI + c e 



1 + d 



+ ... 


(35) 


Using the values (j7]) for A e and A e /, the fit gave b e j = 0.330200(3), c e = 0.7440(4) and 
d = —0.084(3), in agreement with Q and [30]. 


We are now ready to compare our results from the Monte-Carlo simulation with the 
predictions from conformal invariance. We consider the critical Ising model ([3]) in a three 
dimensional cubic lattice excluding all spins outside a sphere of radius R as shown in figure 
[3j The interaction bounds connecting spins inside the sphere with spins outside the sphere 
are also dropped. 

6 We need to be sufficient close to the critical temperature so that the correlation length is much larger 
that the size of the system. 

'The correction proportional to d comes from the leading irrelevant operator in the effective action 
description of the lattice model. 
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Figure 1: Binder cumulant for several system sizes. We used f3 c = 0.22165463 and 1/v = 
3 — A e = 1.58736. The good collapse of the curves shows that (3 C is sufficiently precise for 
the large systems we simulated. 


(e(x)) L 



Figure 2: Expectation values of the energy density operator e(x) at the critical temperature 
of the 3D Ising model with periodic boundary conditions for several linear sizes L of the 
system. The fitting curve is given by the 3-parameter function in equation (35) using the 
values 0 for A e and A e /. 
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Figure 3: Two dimensional cut of a sphere of radius R = 7a embedded in a cubic lattice 
of linear size L = 16. All spins that are outside the sphere are dropped together with their 
interaction bounds (in light grey). 



Figure 4: Expectation value of the energy density operator (minus its value in the infinite 
system) for several values of the radial coordinate r and of the sphere radius R. The collapse 
of all points into a single straight line confirms the prediction (36) of conformal symmetry. 
There are deviations due to finite size effects and due to statistical uncertainty, specially in 
the larger systems. 



















































Figure 5: The same as in figure [4] but with the points grouped into 100 bins with approxi¬ 
mately the same value of r. The statistical error bars are smaller than the size of the dots. 
The black straight line is given by equation (36) with A e from (Jr]) and b ee from © and a t 
from (28). 


In figure [4] we show the one-point function of the energy density operator inside the 
sphere, for various radii R = 7a, 15a, 31a, 63a, 127a. Combining (|6]) with the prediction 
PF1 from conformal invariance, we conclude that 

<■ R /af ' [< e W>B - M = (1 _ + ''' < 36 > 


where the dots stand for terms that vanish in the continuum limit a/R ^ 0 with r/R fixed. 
The plot in figure [5] confirms this prediction and the values of b ee and a e given in [9] and [28] 
In figures [4] and [5] one can notice deviations from spherical symmetry due to the underlying 
cubic lattice, specially for points close to the spherical boundary. For large R we have bigger 
statistical error due to the smaller number of independent samples harvested and because 
the correlation function is multiplied by a large number (i?/a) Ae . This is our first direct 
verification of a non-trivial prediction of conformal invariance. Our second test is related to 
the two-point function inside the sphere. 


From (J5| and (23), we obtain 


(fl/a) 2A " (s(li)s(l 2 )>s = 


(1 - x 2 JR 2 ) A ’( 1 - x 2 /R 2 ) A - 


(37) 


where ( is the conformal invariant ratio introduced in (24) and we neglect terms that vanish 
in the continuum limit. In other words, conformal invariance predicts that the dimensionless 
function 


Fss(xi,x 2 ) 


R 2Aa _ {s(x 1 )s(x 2 ))r _ 

a 2 A a (1 _ x \/R 2 )~ AtT (l — xl/R 2 )~ Aa 


(38) 
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Figure 6: The combination (38) involving the spin-spin two-point function inside a sphere 
against the conformal invariant ratio ( defined in (24). Conformal invariance predicts that 
all points should fall into a single curve up to statistical uncertainties and finite size effects. 


only depends on x\ and x 2 trough the combination In figures [6] and [7J we plot F ss against 
( for many different choices of X\ and x 2 and for several sphere radii. As expected, the points 
collapse in a single smooth curve up to the statistical error bars and finite system size effects. 
Moreover, using 

limF ss (zi,x 2 ) = b 2 srT f aa ( C) (39) 

a —>-0 


and the results (25) and (27) for the asymptotic behaviour of f(Q together with the values 
0.0 and (|28|), we can verify that f aa (Q has the expected asymptotic bahaviour. 


We also preformed a similar analysis for the two point function of the energy density 
operator inside the sphere. In figures [8] and [9j we plot 


p / x _ R 2Ae (e(xi)e( x 2 ))r - (e(x 1 )) R (e(x 2 )) R 

Pee{X 1 ,X 2 )- a2Af ^ _ xlf R 2yA e ^ _ x 2/ R 2yA c 


(40) 


against the conformal invariant (. Using the bulk scaling dimensions [7] and the fact that the 
lowest dimension Z 2 -even surface operator is the displacement operator with Ag = 3, we 
can fit the asymptotic behaviour of the curve to conclude that Cg ~ 0.012. The last value 
is a very crude estimate because F ee has very large statistical uncertainty in the region of 
small (. The main reason for this is that in this region the correlation function is very small 
and it takes a long time to simulate the large systems required to explore the ( —^> 0 limit. 
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Figure 7: The same as in figure [6] but with the points grouped into 100 bins for each system 
size. This reduces the statistical errors and shows better convergence to a single curve. The 
two straight lines are fits to the asymptotic behaviour using the values ([7]), ([9]) and (28). 



Figure 8: The combination (40) involving the connected two-point function of the energy 
density operator inside a sphere against the conformal invariant ratio ( defined in (24). 
Conformal invariance predicts that all points should fall into a single curve up to statistical 
uncertainties and finite size effects. The statistical uncertainty looks biased because we are 
using a logarithmic scale and therefore we can not plot the points with F ee < 0. 
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Figure 9: The same as in figure [8] but with the points grouped into 100 bins for each 
system size. This reduces the statistical errors and shows better convergence to a single 
curve, although the errors are still large for small values of (. The two straight lines are the 
expected asymptotic behaviour using the values (JtJ) , ([9]), (28) and Cq = 0.012. 


5 Conclusion 


We gave strong evidence confirming the non-trivial predictions of conformal symmetry for 
correlation functions of the critical Ising model in a ball geometry. We hope our work 
strengths the confidence in the conformal bootstrap methods that assume conformal sym¬ 
metry from the start. 

It would be nice to obtain more precise measurements of scaling dimensions and OPE 
coefficients of several boundary operators. However, the cubic lattice discretization of the 
ball geometry we used is not ideal for this purpose because it introduces large finite size 
and boundary effects. It would also be interesting to study other BCFT of the critical 
Ising model, like the special and the extraordinary transition. These can be implemented 
introducing another coupling between the boundary spins. Unfortunately, it is not obvious 
how to do this in an elegant fashion in our ball geometry. 

In the absence of boundaries, a fundamental prediction of conformal symmetry is the 
functional form of three-point correlation functions. We plan to verify this prediction with 
Monte-Carlo simulations, in the same spirit of this paper. Such study would also be able 
to check several conformal bootstrap predictions for OPE coefficients of the Ising CFT. It is 
curious that in two dimensions this was done 20 years ago [31]. 
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